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StABUITT  OP  AM  AXISTBBIT  EICAL  C0BFBBSSIB1I  IMfISCID  BAKE 


S.  la.  Gertsenshteyn  and  A.  ¥•  Kasbko 


Sonar  y 

this  article  gives  scic  results  of  the  calculation  of  flow 
stability  in  the  wake  behind  a body  flyiqg  at  a hypersonic  speed.  The 
wave  nuabers,  propagation  rates  and  aaplif icaticn  factors 
characteristic  of  the  oscillations  originating  in  the  wake  are 
obtained.  The  dependences  of  these  values  cn  the  Bach  nusber  and  the 
teopexature  gradient  on  the  axis  and  tie  periphery  of  the  oake  are 
coisidered. 


The  article  also  gives  a short  survey  cf  theoretics'  research  on 
the  stability  of  flou  in  a distant  oaks. 
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1.  Pro*  nany  experinental  studies  [1,  2]  it  is  well-known  that 
sinusoidal  oscillations  originate  in  the  hypersonic  wake  behind  the 
^tkioat"  (at  a distance  cf  10-15  disasters  cf  the  body).  These 
oscillations  develop  rather  quickly,  and  flew  kecoaes  turbulent  at  a 
distance  of  around  25  diaaeters  frea  tic  "throat"  below  the  wake.  The 
study  of  these  oscillations  is  of  special  interest,  and  a nuaber  of 
theoretical  investigations  deal  with  this  question. 

■In  general,  the  results  were  obtained  ior  plane  wakes.  The 
parposc  of  this  investigation  was  to  study  axis yaaetrical 
co spree sib le  wakes. 

Se  will  point  out  that  the  analogy  between  plane  and 
axisyaaetrical  wakes  is  far  fres  ccaplete.  for  exanple,  in  the 
axisyanetrical  case,  it  has  been  shewn  [11  ] that  two-diaeasional 
growing  oscillations  do  not  exist  for  suffioieitly  "saoothed" 
profiles;  only  a three- dine tsicnal  sinusoidal  wave  can  develop  in  the 
flew.  According  to  the  Squire  theorun,  in  the  plane  case,  plane 
disturbances  are  less  stable.  As  the  calculations  show,  the 
corresponding  nuaerical  data  for  axisyaaetrical  and  plane  wakes  also 


' - • 


COC  » 0932  FAGS  3 

differ  considerably  fret  each  ether. 

Obviously,  Helaholtz  was  the  first  to  explain  the  origination  of 
oscillations  in  a wake  ty  the  instability  ci  the  vortex  layers 
feraisg  iq  the  sake.  Rayleigh  [3]  considered  and  solved  the  problea 
of  the  stability  of  a vortex  sheet.  Be  also  solved  the  problea  of  the 
atability  of  a vortex  layer  with  finite  thickness.  The  problea  of  the 
stability  of  flow  in  a wake  with  the  distribution  of  a vortex  with 
finite  intensity  was  first  considered  .by  G.  I.  Petrov  [4].  In  this 
report,  the  velocity  profile  changed  irregularly,  having  five 
branches  (Fig.  1);  the  solution  was  obtained  analogously  to  the  above 
Rayleigh  eolation.  Using  the  relatively  sinple  calculations  in  [4], 
it  was  shown  that,  in  particular,  asynactrical  disturbances  develop 
acre  iateaaively  than  synaetrical.  In  this  sane  report  [4],  the 
generalized  Rayleigh  nethod  was  used  to  consider  the  stability  of  an 
axisyaaetrical  wake  relative  to  zonal  oscillations.  Soae  of  the 
results  of  [4]  were  contradictory,  and  they  were  refined  in  the 
report  by  flichalke  and  Schade  (12).  The  Rayleigh  aethod  was  also  used 
in  £12],  but  the  velocity  ptefile  changed  irregularly  with  nore 
branches  than  in  [4]. 

le  will  point  out  that  the  ovcrwhelning  aajority  of  studies  of 
the  stability  of  wakes  ccacern  plane  wakes,  the  results  on  the 
stability  of  axisyaaetrical  wakes  are  aaialy  given  in  three  studies 
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[11,  4,  13].  As  was  indicated  above,  [4]  ccrsidcrs  the  stability  of 
flow  celative  to  zonal  oscillations.  Scwe  cstisates  and  criteria 
which  are  necessary  and  sufficient  for  the  existence  of  neutral  and 
growing  oscillations  arc  derived  ip  [11]  analogously  to  the  known 
criteria  for  the  plane  case.  The  necessary  condition  for  the 
existence  of  growing  disturbances  cbtained  in  [11]  is  the  nost 
interesting.  Below  we  will  briefly  explain  bow  this  condition  can  be 
obtained  and  what  ccnclasicrs  can  easily  be  drain  fron  it. 


1JIW" 


,!F'"T 


The  problen  of  the  stability  cf  ap  axisyanetrical  wake  in  a 
perfect  ipcoapressible  fluid  can  be  reduced  to  the  problen  of  the 
proper  values  for  the  scccnd-order  eguatiop 


Bare  the  boundary  conditions  are  the  conditions  of  the  linitedness  of 
the  solution  at  « - o and  «— «•  . 0 (o  is  the  velocity  profile  of 

the  nain  flow,  c and  n are  the  assigned  percenters  (the  disturbance 

was  represented  as: 

$ *•  g (*)  **?(**  * * ♦ »*  f)  . 

Multiplying  eguation  (1)  by  ~ (here  G is  the  fnnetion  conjugate 

w*  C 

to  G)*  integrating  the  expression  cbtained  irca  v-o  to  (-•  , 
integrating  by  parts  using  the  known  boundary  conditions,  and 
separating  the  real  and  inaginary  part,  it  is  easy  to  obtain: 
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Wi  • o 


Here 

*W-Frr*Wi  e--^  • 

If  C-  * 0#  then  it  follcws  froa  (2)  that  Q*  shcald  change  signs  in 
the  interval  (0,  •) , i.c.f  at  C;  1 0,  C'  sfcculd  return  to  zero  at  a 
certain  point  in  interval  ((,  •).  Frca  this  necessary  condition,  it 
is  not  hard  to  show  that  C's  Oat  n * 0 fci  o(t)*o,/(i*^)1  , 

for  enaaple. 

Veport  [11]  shows  that  for  this  prefile,  growing  disturbances 
can  only  exist  at  n - 1,  i.e.,  only  "sinusoidal*1  oscillations  can  be 
unstable. 

Article  [11]  points  out  that  this  cooclusicn  can  only  be  used 
for  "saoothed”  velocity  pre files. 

For  profiles  with  a narked  change  in  velocity,  this  conclusion 
is  knpwa  to  be  incorrect.  lie  specific  conpctations  for  growing 
disturbances  are  not  given  in  [11]. 

Study  [13]  aainly  ccnccrns  plane  vakea  in  a coapressible 
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inviscid  fluid.  For  the  axisyaaetrical  ca$f,  this  report  gives  only 
iKlattd  calculations  fee  the  Causa ian  distribution  of  the  velocity 
profile  V(t j.f**'  at  a * 1,  i.e.,  for  "sinusoidal"  waves.  Report  [13] 
not  oaly  gives  specific  calculations,  hot  certain  necessary  criteria 
for  the  existence  of  neutral  disturbances  aid  a certain  sufficient 
condition  for  the  existence  of  growing  disturbances  are  derived.  This 
study  is  analogous  to  the  studies  conducted  for  an  incaa pressible 
fluid  [5,  6].  The  authors  also  used  the  studies  of  Lin*  and  L.  Lix  of 
the  *tability  of  a coapressible  boundary  layer,  which  have  not  been 
published  in  the  Soviet  Onicn,  either. 


Jcsides  the  above  study  [13],  the  stability  of  a coapressible 
wake  sas  also  studied  in  [7,  It,  15].  The  plane  case  is  investigated 
ia  [7,  It].  Article  [15]  ccrsiders  the  plate  and  axisyaaetrical 
cases*.  Bare  the  sain  flew  was  assuxed  to  be  ccrstant  in  the  vicinity 
of  the  axis,  i.e.,  the  stability  cf  an  axisyaaetrical  vortex  sheet  in 
a coapressible  gas  was  considered. 


2.  It  is  convenient  to  consider  the  initial  systea  of  equations 
it  the  cylindrical  coordinate  systea  [8]  fer  studying  oscillations  in 

aa  axisyaaetrical  sake: 


T 


H *»*~S7  *,4~St  '7 

I »r,  » *v  / ** 

i-nvi?V7jF*'jr'«#3r--*  * 

^ , ,1^*,  .fT; , (3) 

J*  » f • i»  » Jt  f 

*1  * *L5?  » *«K ) . .&.o 

it  it  » 39  it  * 

4T  P df 

c*>-iVk7Tt 

P • fpr. 

We  vill  represent  the  unknown  soluticn  cf  the  systea  in  the  fora  of 
the  sun  of  the  sain  fled 


V (■»  « V » • 0 v * 0 

"«», ' l»  l v*(t)  • •*»■ « * *»•*  • 

r ■ f ^ f p • S ) 

aad  a certain  sufficiently  snail  disturbance 


• . • V ♦ »'  , r • w'  v • v' 

• ■ a • * t • • * • 


P-P  * p' 


r . f ♦ r'  . 


Substituting  (4)  is  (3)  and  disregarding  the  teras  which  are 
guadratic  relative  to  the  disturbance  ccaycnents,  we  will  obtain  a 


i 


: 


system  of  six  linear  equations  relative  to  six  unknown  functions 
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4* 

V ^ * 

<>*• 

■S 

• Jl 

4 L i 

t t 

4< 

P 'f  'T  * f T’  . 

If  ve  search  for  all  of  tie  disturtance  conjonents  in  the  forn 

J * ^ (*)  * >atl  * iff]  , 

it  is  easy  to  obtain  the  systen  of  ordinary  differential  equations 
for  determining  function  j (x)  frcn  systen  ff)  : 
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* * r + • 

i*  *,  * vt *,  .-I 

-»  ,,  * i . * 

</»  V,  ♦V,  t W «>,  « - — i P 

*? 


W*!**!  Vt  (*)*i<*Vt(t)V‘l  •- J- i*  P (< 

A 

• . r* ' t . + * v.  7 

♦ *- v;  j>  ♦ -*/»;  ♦*<**’■,  *T  \ *o 

pifiT  *p»tT'(x)*pi«VtT  =(t-rt)p[i'trl*  »t.il  J 

P-  f?  * ?T  . 


*e  will  introduce  the  diaensic nless  variables: 


f - 


(7) 


In  diaensionless  variables,  system  (6) 


V*  «%♦!', 


tj>?,  *Vf  i<* 


A 

r it*P 

5 fl'r 


has  the  fora 

(a) 

(0) 

do) 
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v*p  * 


iV  • 


(il) 


'('-r.)plu*t*K  •**  *•  ] 


(12) 


(13) 


• n 

Bexe  B is  the  Hach  nuaber:  « . 

Pox  farther  calculations,  it  is  cotvcnieqt  tc  reduce  systea  (8)  to 

A. 

one  equation.  He  will  write  the  equation  fex  pressure  P.  First  it  is 

* A A 

expedient  to  find  the  expression  fer  »,  fica  E and  P*.  In  order  to 
do  this,  it  suffices  to  substitute  the  expression  for  v,  froa 

A 

equation  (8)  in  equation  (1C).  Substituting  the  expression  for  T froa 
equatioa  (13)  in  equation  (12),  then  expressing  it  by  and 

obtained  froa  (11)  instead  of  £ we  will  ettain  the  relationship 
aacag  P e, . . If  we  now  ass  the  previously 

deter ained  expressions  for  »,  ,e#  and  »,  by  t ead  P* , we  will  obtain 

A 

the  following  equation  for  E: 


(14) 


mmmmm 


-jr-  tmm 
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Here 


3.  It  is  easy  to  substitute  tbs  boundary  conditions  for  equation 
(14):  1(0)  * 0 at  *-o,  ud  at  «->•  , the  absolute  value  of 
solution  *(*)  should  be  liaited.  However,  it  is  not  advisable  to 
aunerically  integrate  equation  (14)  near  pcint  * -0,  since 
expressions  cf  the  fora  *n<3  T>  *re  in  tic  coefficients. 

Therefore,  it  is  logical  to  use  the  asyaptctic  representation  of  the 
soluticn  in  the  vicinity  cf  zero. 


In  order  to  perfora  the  nuaericel  calculation,  it  is  also 
necessary  to  replace  the  boundary  ccnditicy  at  infinity  by  the 
boundary  condition  at  a finite  distance.  This  can  be  done  by  using 
hycun  asyaptotic  expansions.  He  will  give  tie  corresponding 
computations  below. 


If  v,-a  i and  fa  i at  , equation  (14)  ass  uses  the  foru: 


«»  i 


.<>-)  M*  'if/**)  P T?  * 


0 . 


(18) 


cr 
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**•  asyaptotic  behavior  of  the  Halted  eoletico  to  equation  (16)  at 
infinity  la  nail-known  [9,  10]; 


P(%) 5 * 

n 


where 


A*  t * o 


le  are  not  considering  the  case  Hex  * 0,  since  the  boundary  condition 

at  infinity  is  satisfied  autoaatically  for  any  £ at 

*Yrf  ♦//-•<*«  o (the  case  of  a continuous  spectrua)  . Using 

^ to  see  that  at  sufficiently  large  values  of  r,  the 

following  relationship  should  be  satisfied  kith  great  accuracy: 


. (it) 

relationship  (19)  will  also  replace  the  boundary  condition  at 
infinity  in  the  future. 
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1«  will  Kit*  the  asyaptotic  solution  to  aquation  (It)  near 
teto.  It  is  convenient  tc  rewrite  equation  <m)  as  follows: 


(30) 


■e  will  represent  equation  (20)  in  the  fora  of  a series: 

Gather  lag  the  terns  of  the  fora  *w,  «e  will  have: 

\(\.  t)ca°  X(-Ca)*  j'ca 

She  nee: 

*•  * 

* 7 • 

It  is  easp  to  find  the  reaaining  cccfficients  cf  the  series  froa  C0. 
For  epaaple,  by  gathering  the  teras  of  the  lorn  «*"  , we  will  obtain 


C, ■* «f I 1 

rmt  ■ #te#  J • 


thus,  the  uahaoaa  solution  behaves  like  «r  near  sero.  Here  it 
is  assuaed  that  r+  o . For  f-o  , we  can  slec  obtain  the 
corresponding  asyaptotic  representation  analogously  to  the  above 
procedure  [27]. 
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4.  These  calculations  Bade  it  possible  to  get  an  idea  of  the 
nature  of  the  oscillations  in  a supersonic  take.  As  an  exaaple,  ue 
considered  flow  in  the  wake  tebind  a tody  it  tie  range  of  Rach 
nuaters  fros  10  to  30.  Ihe  flew  instability  is  considered  in  the  wake 
behind  the  "throat,"  at  a distance  cf  10-15  tines  the  size  of  the 
body.  Ihe  data  on  the  velocity  and  teaperatere  distributions  were 
taken  iron  [ 1,  2,  13}: 


■ere  u(t)  is  the  diseasicnless  velocity,  T,\>  is  the  teaperature 
distribution,  *,(•)  is  the  value  of  the  enthalpy  in  the  cross  section 

we  selected  on  tbe  axis  of  the  wake,  e*d  e„  is  the  value  of  the 
enthalpy  in  tbe  inccsing  flew. 


The  results  of  the  calculaticca  are  partially  shown  in  graphs 
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Da  rigoras  1,  2,  3 and  4 skov,  tka  aaplif icatioa  factor 
dacxaaaaa  aitk  tka  incraaaa  in  tka  Hack  nualer  fl,  bat  extrenely 
slightly.  Therefore,  tka  aall-kacva  aarkad  incraaaa  in  flow  stability 
aitk  tka  incraaaa  in  tka  Hack  nnabcr  can  neat  probably  ba  explainad 
by  tka  kigk  phase  velocity  cf  tke  propagaticn  of  the  disturbance  in 
tka  ooapressible  wake  (fig.  5)  - (the  distar  tare as  are  polled  into 
tka  flow  vary  quickly*  vitkout  being  alia  tc  develop). 


‘H.!.  I jr  -W'jm* 
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rig. 


5. 


Jt  is  interesting  to  ccapare  the  resalts  obtained  with  the  known 
remit*  [13]  of  the  stability  of  flat  jets*  1h«  comparison  shows  that 
the  aaplif ication  factor*  for  flat  jets  are  approxiaately  four  tises 
higher  than  for  ax isyanetrical  jet*. 

Vroa  these  calculation*  we  can  also  see  that  the  extreae  wave 
naaber  •*  , which  corresponds  to  the  naxiau*  aaplif icaticq  factor, 
decreases  with  the  incxeaic  in  the  Bach  nuaber,  while  the  phase 
velocity  cf*  reaains  virtually  unchanged.  I*  will  also  point  oat  that 
the  phase  velocity  varies  little  with  the  change  in  the  save  nuaber  * 
(at  a fixed  Bach  nuaber)  - fig.  5. 
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Separate  calculations  fcr  confuting  the  effects  of  a hot  "core" 
fct  aa  axiaynnetcical  wake  were  Bade  in  this  study. 

•he  aaplif ication  factcr  «c*  increases  approxiaately  1.5  tiaes 
(Pig.  6)  when  AT  increased  (froa  0.1  tc  2.0),  while  the  extrene  wave 
cuaber  a*  decreased  by  appioxiaately  1/10  cf  its  initial  value  when 
AT  increased  froa  0.1  tc  2.C.  The  rate  cf  propagation  of  the 
disturbance  in  the  range  cf  change  ia  question  changes  very 
insignificantly.  He  will  pcint  out  that  for  observers  aoving  at  the 
rate  of  propagation  of  the  disturbance  , the  local  velocity 

cf  the  flow  whose  stability  is  being  investigated  is  lower  than  the 
local  speed  cf  sound.  This  situaticn  hclds  throughout  the  range  of 
change  in  AT  and  H in  guesticn. 
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